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Abstract 

Measurements of charged pion and kaon production are presented in centrality selected Pb+Pb 
collisions at AOA GeV and 158A GeV beam energy as well as in semi-central C+C and Si+Si 
interactions at AOA GeV. Transverse mass spectra, rapidity spectra and total yields are determined 
as a function of centrality. The system-size and centrality dependence of relative strangeness 
production in nucleus-nucleus collisions at 40^4 GeV and 158^4 GeV beam energy are derived from 
the data presented here and published data for C-l-C and Si-|-Si collisions at 158^ GeV beam 
energy. At both energies a steep increase with centrality is observed for small systems followed by 
a weak rise or even saturation for higher centralities. This behavior is compared to calculations 
using transport models (UrQMD and HSD), a percolation model and the core-corona approach. 

PACS numbers: 25.75.-q 
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I. INTRODUCTION 



Strangeness production has always been an important observable for probing the state of 
matter created in heavy-ion collisions. The yields of the various strange hadrons, properly 
normalized to those in elementary nucleon-nucleon interactions, are expected to exhibit dif- 
ferent behavior in partonic than in hadronic environments The number of created strange 
quark-antiquark pairs per participating nucleon is significantly higher in central Pb+Pb col- 
lisions than in p+p interactions at all collision energies. A maximum in the yield of the 
K'^ relative to that of the tt"'" in central Pb+Pb collisions was observed by NA49 at about 
30A GeV beam energy and interpreted as evidence for the onset of deconfinement. 

Microscopic models such as HSD and UrQMD do not predict a maximum. In the statisti- 
cal model this structure can only be reproduced by including additional features like those 
introduced in j4|, i.e. high mass resonances beyond the ones listed in jsl and a significant 
contribution by a meson states. There are no indications for a similar non-monotonic energy 
dependence of this ratio in p+p interactions. In order to better understand the observed 
maximum in the energy dependence observed in central Pb+Pb collisions, measurements of 
the system size dependence of strangeness production are needed for different strangeness 
carriers and at various energies. 

NA49 has published data on the system-size dependence of hyperon (A and S) 6| and 
(anti-)proton yields and distributions 0] in Pb+Pb collisions at 40^4 GeV and 158^4 GeV 
beam energy. Results on pions, kaons, mesons and A hyperons in C+C and Si+Si at 



158y4 GeV were reported in Ref. [8[. Pion and kaon production was also studied in p+p 
interactions at 158 GeV/c [q, [io|. The present study completes these results by a mea- 
surement of the centrality dependence of pion and kaon production in Pb+Pb collisions at 
AOA GeV and 158A GeV beam energy as well as for C+C (66% most central) and Si+Si 
(29% most central) interactions at 40^4 GeV. These results provide a link between the cen- 
trality selected data on pion and kaon production at lower (SIS, AGS) and higher (RHIC) 
energy. At SIS and AGS energies kaon and pion production was studied in centrality se- 



lected Au+Au collisions as well as in collisions of smaller systems such as Si+Al |lll-ll3|. 
Here, a rather linear increase is observed for relative kaon production with the number of 
participants. On the other hand, at RHIC energies a clear saturation behavior is found for 



collisions with more than approximately 60 nucleons participating in the collision 
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The same observation has a. 
Si+Si and Pb+Pb coUisions 



Statistical models based on the grand-canonical ensemble (G-CE) have been used to de- 
scribe strange particle yields in central collisions of heavy nuclei at high energies 
For smaller reaction volumes, because of the lower number of produced strange quarks, 
restrictions due to strangeness conservation arise, which suppress strangeness production 
with respect to the infinite volume limit. This is taken into account by using canoni- 
cal ensembles in the statistical model calculations. The transition from the canonical to 
the grand-canonical ensemble would thus result in a characteristic volume dependence of 
strangeness production 18|, ll9|]. When comparing statistical model results to experimental 
data in terms of reaction volume, the latter has to be related to an experimental observable. 
In Refs. 19|, |20| the reaction volume was assumed to be proportional to the number of 



1 this approach the data at SPS energies 



21| . Good agreement can be achieved, if 



wounded nucleons (see section III). However wit 
and above could not be described satisfactorily 
the relevant volume is subject to additional geometric constraints which restrict strangeness 
conservation to subvolumes, as realized for example in the percolation models. The corre- 
sponding calculations suggest that not only for small systems several subvolumes are formed, 
but that also in central collisions not all nucleon-nucleon collisions are included in the central 
volume 
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model 



A simpler implementation of this ansatz is realized in the so-called core-corona 



23| . It is built on the assumption that the fireball created in nuclear collisions 



is composed of a core, which has the same properties as a very central collision system, 
and a corona, which is a superposition of independent nucleon-nucleon interactions. The 
sizes of core and corona can be defined by the condition that the corona is formed by those 
nucleon-nucleon collisions in which both partners interact only once during the whole col- 
lision process. The fraction of single scatterings is calculated using straight line geometry 



n 

as described in the Glauber model [24]. This approach was recently success 



■ul 



y applied to 



the system-Size depeudeuce of strangeness production at RHIC energies and for the 

above mentioned hyperon production at SPS [6J]. The common features of such model calcu- 
lations at fixed energy are a fast increase of relative strangeness production with system size 
for small reaction volumes (below approximately 60 participating nucleons) and eventual 
saturation for large system sizes. 



II. THE NA49 EXPERIMENT 



The NA49 detector is a large acceptance hadron spectrometer at the CERN SPS 27|] . The 
main components are four large time projection chambers (TPCs) and two super-conducting 
dipole magnets with a one meter vertical gap, aligned in a row, and a total bending power 
of 9 Tm. The magnetic field was set a factor of four lower for data taking at 40^4 GeV 
than at 158A GeV. Two two-meter long TPCs (VTPCs) inside the magnets each with 72 
pad-rows along the beam direction allow for precise tracking, momentum determination, 
vertex reconstruction, and particle identification (PID) by the measurement of the energy 
loss (dE/dx) in the detector gas with a resolution of 6%. The other two TPCs (MTPCs) 
downstream of the magnets have large dimensions (4 m x 4 m x 1.2 m, 90 pad-rows) and 
provide additional momentum resolution for high momentum particles as well as PID by 
dE/dx measurement with a resolution of around 4%. A momentum resolution in the range 
cr{p)/p'^ = (0.3 — 7) ■ 10^^ (GeV/c)~^ is achieved. Two time-of-flight scintillator arrays of 
891 pixels each (TOF), situated just behind the MTPCs symmetrically on either side of 
the beam axis, add additional K/vr separation power in the laboratory momentum range 
from one to ten GeV/c (near mid-rapidity for kaons). A veto-calorimeter (VCAL) placed 
20 meters downstream of the target accepts all beam particles and projectile fragments 



as well as most of the spectator neutrons and protons 27|]. The geometrical acceptance of 



VCAL was readjusted for each beam energy by means of a collimator for optimum projectile 



spectator coverage. VCAL is used for off-line (on-line) selection of event centrality in P 



b+Pb 



(C+C and Si-|-Si) collisions. The NA49 detector is described in detail in reference [27 1. 

For Pb+Pb collisions a primary SPS beam of 40^4 GeV and 158^4 GeV was directed onto 
a Pb target while for C+C and Si+Si interactions a fragmented Pb beam of 40A GeV was 
used. The fragments were identified by magnetic rigidity {Z/A = 0.5) and by a pulse-height 
measurement in a scintillation counter in the beam. For C+C and Si+Si interactions beam 
particles with Z around 6 and Z around 14, respectively, were selected online. A refined 
Z selection was done offline, which yielded a pure C-beam, whereas for the "Si-beam" the 
different charges could not be separated well enough, thus a mixture of ions with Z = 13, 
14, and 15 (intensity ratio 34:44:22) were accepted. The target disks had thicknesses of 
10 mm (carbon), 5 mm (silicon) and 200 /xm (lead) corresponding to 7.9%, 4.4%, and 0.5% 
interaction lengths for the respective beam particles. 



The Pb (light ion) beam passed through a quartz Cherenkov (scintillator) detector, from 
which the start signal for the time-of-flight measurement was obtained, and three stations of 
multi-wire proportional chambers which measured the trajectory and energy loss of individ- 
ual beam particles. A minimum bias trigger for Pb+Pb interactions was derived from the 
signal of a gas Cherenkov device right behind the target. Only interactions which reduce the 
beam charge and thus the signal seen by this detector by at least 10% were accepted. The 
interaction cross section thus defined is 5.7 b at both energies. The trigger for the C-l-C and 
Si+Si data taking at AOA GeV was based on the energy deposited in VCAL. The 66% most 
central C+C and 29% most central Si+Si interactions were selected by requiring an energy 
in VCAL below the imposed threshold. These fractions are determined using the inelastic 
A+A cross sections (10% uncertainty), the trigger cross sections (5% uncertainty), and a 
simulation of the VCAL response. 



III. DATA ANALYSIS 

The data sample used for this analysis was collected in specific minimum bias runs of 
the NA49 experiment and thus are different from the central data presented earlier [2]. The 
recorded minimum bias Pb+Pb collisions were divided into five consecutive centrality bins: 
CO - C4 (see Tab. [land Ref. (28|). The centrahty selection is based on the forward going 
energy of projectile spectators as measured in VCAL. We quantify centrality by the fraction 
of cross-section according to intervals of forward going spectator energy. For each centrality 
interval a characteristic quantity, the mean value of "wounded nucleons" (A^w) is calculated. 
A nucleon is considered wounded, if its interaction occurs in the nuclear overlap volume. 
(Aw) is determined by generating VCAL spectra with events from the VENUS 4.12 Monte 



Carlo code 29|. The simulation took into account the energy resolution of the calorimeter 
and contributions of participants to the energy recorded by VCAL. After cross calibration 
of experimental and simulated spectra for effects of the experimental trigger in the most 
peripheral centrality bin the meqan number of interacting (wounded) nucleons ((N^)) was 
extracted from the model data for each of the selected cross section fractions [30|. The 
latter were selected to be identical at 40^4 GeV and 158^4 GeV. Thus only insignificant 
differences of {N„) are observed between both energies except for the 10% discrepancy in 
the most peripheral centrality interval which we attribute to the slightly different online 
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C2 
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212 ±3 


34 


40A 


C3 
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144 ±4 
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33.5-43.5 


93 ±7 
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0-5 


357 ±1 
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158^ 


CI 


5-12.5 


288 ±2 


24 


158A 


C2 


12.5-23.5 


211 ±3 


37 


158A 


C3 


23.5-33.5 


146 ±4 


35 


158A 


C4 


33.5-43.5 


85 ±7 


35 



TABLE I. Overview of the analyzed data and centrality classes for minimum bias Pb±Pb collisions. 
The centrality is given as fraction of the total inelastic cross section. (A^w) is the average number 
of wounded nucleons per event, and Nevent the number of analyzed events. 

trigger conditions. Table H] summarizes these numbers together with the number of analyzed 
events. For convenience we quote here the centrality of the C±C and Si±Si collisions at 158 
A GeV ^ which are used for comparison in this paper: the most central 15.3% ((A^w) = 
14±2) and 12.2% ((A^^) = 37±3), respectively. 

The recorded Pb±Pb events are contaminated by background due to non-target inter- 
actions. This contamination was minimized by cuts on the fitted vertex position and the 
quality of the fit. It is negligible for near central collisions and amounts to less than 5% for 
the most peripheral collisions (see Similar cuts were applied to the C±C and Si±Si 



collision events. Details for the latter are given in Table 4.1 of 31|. 

The track finding efficiency and dE/dx resolution were optimized by track quality criteria, 
some of which are different for pions and kaons. These criteria are meant to select phase space 
regions with low track densities, where the background and efciency corrections are small 



and approximately uniform. Only those particles were considered for which the transverse 
momentum kick due to the magnetic field points into the same direction as the azimuthal 
angle projected onto the deflection plane. In addition the number of reconstructed points 
on the particle track had to exceed 50% of all measurable points. Further track cuts on 
azimuthal angles with respect to the bending plane and for a minimum of measured and/or 
measurable points are meant to obtain high quality track samples. These track cuts were 
chosen differently for different particle species and data samples to optimize phase space 
coverage: For negatively charged hadrons and kaons azimuthal angles within ±30 degrees 
with respect to the bending plane were selected except for kaons in the C+C and Si+Si 
events where ±45 degrees were allowed to enhance statistics. Tracks of negatively charged 
hadrons in C±C and Si±Si events were required to have at least 10 potential points in the 
TPCs, while in Pb±Pb events a minimum of 30 measured points was asked for. Kaon tracks 
in C±C and Si±Si events had to have a minimum of 50 measured points in all TPCs and 
at least 1 potential point in the VTPCs. In Pb±Pb events they were required to have a 
minimum of 50 measured points in the MTPCs and at least 1 measured point in the VTPCs, 
if the number of measurable points was larger than 10. Finally kaon candidate tracks in 
the sample of Pb±Pb events were required to point to the main interaction vertex with a 
precision of 2 cm in the x- (bending) and 1 cm in the y-direction. 

While kaon spectra were determined by means of particle identification, pion spectra 
were obtained by subtracting from all negatively charged particles the contributions from 
electrons, kaons, and anti-protons (see below). Kaons were identified using the energy loss 
[dE/dx) information from the TPCs and time-of- flight measurements in the phase space 
covered by the TOF walls. For the analysis based on dE/dx information alone, raw yields 
of K"*" and K~ were extracted from fits to the dE/dx distributions in bins of laboratory (piab) 
momentum (11 bins from 5 to 50 GeV/c in logarithmic scale for the 158^4 GeV data set and 
10 bins from 4 to 32 GeV/c in logarithmic scale for the 40^4 GeV data set), and of trans- 
verse momentum (px) in 20 bins from to 2 GeV/c in linear scale. See |32| for a detailed 
description of the fitting method. The energy loss information from the MTPCs (VTPCs 
and MTPCs) was used in the analysis of minimum bias Pb±Pb (C±C and Si±Si) collisions. 
Bins with less than 500 particles were not used in the further analysis to ensure sufficient 
statistics for the fits to the dE/dx distributions. This cut limits the accessible range in 
Pt to values below 1 to 1.5 GeV/c depending on piab- The phase space in piab and thus the 
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acceptance in rapidity y is limited by the fact that only for momenta from piab = 4 GeV/c to 
Piab = 50 GeV/c the mean energy loss of kaons is sufficiently separated from those of pions 
and protons to allow for the extraction of their yield by unfolding of the dE/dx distributions. 
This restriction limits the acceptance of kaons to forward rapidities. The resulting distribu- 
tions were corrected for geometrical acceptance, in-flight decay of kaons, and reconstruction 
efficiency. The flrst two corrections were determined by detector simulations using GEANT 



33| . The reconstruction efficiency is track density dependent only for Pb+Pb collisions. It 



was shown earlier [8| that the corrections are negligible for C+C and Si+Si interactions. 
The correction factor was determined by embedding simulated tracks into real events in 
three centrality bins ranging from (0-12.5) %, (12.5-33.5) %, and > 33.5 %. Track losses are 
of order (5-10) % for central collisions at 158^4 GeV and drop to only (1-2) % in the most 
peripheral collisions and for 40^4 GeV beam energy. Systematic errors were determined by 
varying track selection criteria and the parameters in the dE/dx unfolding procedure. In 
particular we decreased the minimum number of measured points from 50 to 30 and used an 
azimuthal cut of ±50 instead of ±30 degrees. The asymmetry parameter in the dE/dx un- 
folding procedure, which accounts for the asymmetry of the distributions of the truncated 
mean (for details see 32|]), was varied within the limits 0.6 and 0.8. Different combinations 
of these parameter settings resulted in the systematic errors indicated in Tables II and III. 
The spectra at mid-rapidity were obtained using the combined dE/dx and TOF informa- 
tion (Pb±Pb only; statistics were not sufficient for C±C and Si±Si). The distributions of 
TV? calculated from the recorded time of flight (cttof ~ 65 ps) were parametrized in bins 
of laboratory momentum as the sum of three Gaussians for pions, kaons and protons plus 
exponential tails stemming from unresolved double hits or coincidence with 7 conversion in 
the scinitillators. For laboratory momenta below 2.5 GeV/c, this measurement is sufficient 
to separate kaons from pions and protons with high efficiency (~ 95%) and high purity 
(> 96%). For higher momenta, the independent measurements of time of flight and specific 
energy loss are used simultaneously. In the two-dimensional plane of and dE/dx, kaons 
were identified by an elliptical region centred at the expectation values and with half-axes 
equal to 2.5 a^'^ and 2.5 o"dE/dx; respectively. In addition, regions where the local kaon con- 
tribution, according to the parametrisations of m? and dE/dx, is below 0.7 were excluded. 
Each track in the kaon identification area entered the transverse monentum spectrum with a 
specific weight accounting for the efficiency of particle identification and for the efficiency of 
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the time-of-flight detector. The identification efficiency is defined by the ratio of the integral 
of the two-dimensional kaon distribution in and dE/dx inside the identification area to 
its total integral. It does not depend on the centrality class and is determined as function 
of track momentum. For the kaon identification employed here, the efficiency is about 95% 
and varies only slightly with momentum. The contamination by pions and protons is below 
1% for all momenta and is not corrected for. Efficiency losses of the time-of-fiight detectors 
are caused by multiple hits in one scintillator and quality cuts on the deposited charge. 
They were determined from the data themselves for each channel and in each centrality 
class separately. The average efficiency is 82% and varies slightly with centrality. Finally, 
the transverse mass spectra were corrected for the geometrical acceptance and the in-fiight 
kaon decay, which was obtained by a detector simulation with GEANT. 

The separation of negatively charged pions from kaons and anti-protons in the acceptance 
of the NA49 detectors by means of the energy loss measurement is limited to high transverse 
momenta (px larger than 500 MeV/c at mid-rapidity and larger than 200 MeV/c at ?/cm=l 
for 158^4 GeV beam energy). Therefore a different method was chosen for the evaluation 
of vr" production (see jsl). Yields of all negatively charged hadrons (h~ from the primary 
vertex were extracted in bins of rapidity and assuming the pion mass. The contamination 
by K~, p, e~, and by particles from the decay of strange particles and from secondary 
interactions, both reconstructed at the primary vertex, was estimated by simulation and 
subtracted. The calculation of theses corrections is based on events generated by the VENUS 
model [29|, propagated through the NA49 detector using GEANT, and processed by the 
standard reconstruction chain. The ratio of reconstructed to all simulated particles yields 
correction factors in all bins of rapidity and which account for all background tracks in 
the negatively charged hadron yields. Care was taken to account for differences between 
simulated and measured hadron yields by appropriate scaling factors. These were obtained 
from a comparison of the simulated yields with measured data. Where measured yields 
are not available, particle ratios from statistical model calculations were used to infer the 
corresponding yields. In general distributions of background particle from VENUS events 
were scaled by a global factor. The largest corrections to the pion yields arise from K~ 
in Pb+Pb collisions (158^4 GeV) at high px (30% for px > 1-2 GeV/c) where pion yields 
are low and kaons prevail. The K~ distributions were therefore scaled differentially (and 
iteratively) in rapidity and px in order to reproduce the shapes of the measured spectra. The 
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differences in the correction factors on the negatively charged particles from K~ was found 
to be in the range of 1-3%. As this is a very small change, we found it justified to use global 
factors for the other contributions. At mid-rapidity and low (< 0.1 GeV/c) corrections 
amount to 30% - 40% due to contributions from (secondary) electrons (158^4 GeV). This 
large correction requires a careful evaluation as the contribution relies fully on the GEANT 
calculation. We therefore identified and removed electrons from the track samples in both 
the simulation and in the data (for the latter based on dE/dx measurements). Results agree 
within 5% - 10% in the corresponding bins of rapidity and transverse momentum and within 
2% - 3% for the overall pion yield. For the overall evaluation of systematic errors on the pion 
yield we performed further analyses with different track selection parameters (minimum of 
50 instead of 30 measured points), reversed magnetic field, and a different data set with 
slightly varied settings for the online selection of minimum bias events. Variations of the 
pion yield are less than 5% in centrality bins C0-C3 and 10% in bin C4. This higher error 
in the latter bin is due to the uncertainty of the event selection for the most peripheral bin. 
In C+C and Si+Si collisions at 40^4 GeV vr" were analyzed with both methods (h~ and 
dE/dx). The differences were found to be below 5%. 

The subtraction method is not well suited for the determination of the 7r+ yields at our en- 
ergies, because positively charged particles have a large contribution of protons. However, for 
Pb+Pb collisions yields at mid-rapidity and full phase space can be inferred using the feed- 
down corrected tt^/ti^ ratio from the combined dE/dx and TOF analysis in an acceptance 
which is a banana shaped region in the rapidity-px plane within < J9t < 1400 MeV/c and 
2.1 <y < 3.1 (0 < pt < 1300 MeV/c and 2.5 < y < 5.0) at AOA GeV (158A GeV). This 
method relies on the assumption that the pt and rapidity distributions of vr" and vr"*" have 
the same shapes. In VENUS events this assumption leads to differences of 2% between the 
direct and scaled 7r+ multiplicities. The extracted /n^ ratios are 0.93 ± 0.2 and 0.90 ± 0.1 
at 158A GeV and iOA GeV beam energies, respectively and are subject to a 4% systematic 
uncertainty, which is derived from the difference between two different methods (see and 
j^). We assume that this ratio is centrality independent. For C+C and Si+Si collisions 
the 7r~/7r+ ratios were determined at 158^4 GeV beam energy (0.99 and 1.02) in [8j]. We 
use these same values at 40^4 GeV. More details on the analysis of minimum bias Pb+Pb 

. The C+C and Si+Si analyses of pion and kaon yields at 
and 



collisions can be found in 
40^4 GeV are described in 



35[, respectively. 
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IV. RESULTS 



A. Transverse mass spectra 
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FIG. 1. (Color online) Transverse mass spectra at mid-rapidity from the combined dE/dx and TOF 
analysis for K+ and in Pb+Pb collisions at AOA GeV (-0.1 <y < 0.1) and 158A GeV (-0.2 < 
y < 0.2) beam energy. Lines show exponential fits for 0.0 GeV/c^ < uit — tuq < 0.8 GeV/c^ 
(40^ GeV) and 0.0 GeV/c^ < mx - mo < 1.0 GeV/c^ {158A GeV). Different centralities are 
scaled for clarity. The error bars indicate the statistical uncertainty. 




FIG. 2. (Color online) Transverse mass spectra at mid-rapidity (0 < y < 0.2) for vr^ in Pb-|-Pb 
collisions at AOA GeV and 158^4 GeV beam energy. Lines show exponential fits for 0.2 GeV/c^ 
< niT — rno < 0.7 GeV/c^. Different centralities are scaled for clarity. The error bars indicate the 
statistical uncertainty. 
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FIG. 3. (Color online) Transverse mass spectra of tt^ (0 < y < 0.2) and kaons (0.8 < y < 1.0) at 
40^4 GeV in C+C and Si+Si collisions. Lines show exponential fits for 0.2 GeV/c^ < mx — niQ < 
0.8 GeV/c^ (■^~^) 3.nd 0.1 GeV/c^ < mx — mo < 0.5 GeV/c^ (kaons). The inverse slope parameters 
for K+(K") are 171 ±18 MeV for C+C and 203 ±10 MeV for Si+Si collisions (137 ±24 MeV for 
C+C and 166 ±13 MeV for Si+Si collisions). The inverse slopes of pions are given in Table HVl 
The error bars indicate the statistical uncertainty. 

Transverse mass spectra at mid-rapidity are shown in Figs. [1] and [2] for kaons (combined 
dE/dx and TOF analysis) and pions for different centrality classes in Pb+Pb collisions at 
AO A GeV and 158 A GeV beam energy. Fig. [3] presents results for kaons and pions near mid- 
rapidity in C+C and Si+Si collisions at AOA GeV. The mx spectra (mx = a/ (m^ + p|)) 
were fitted by an exponential function with inverse slope T 



dN 



C ■ exp ^- 



' T ) 



(1) 



mTdrriT dy 

in mx — rriQ ranges as given in the figure captions. In Pb+Pb collisions the transverse mass 
spectra of kaons are well described by this functional form, while pion spectra deviate at 
high and low transverse masses. The mx spectra of kaons in the light systems exhibit a two 
to three times statistical error downward deviation at low mx. For a model independent 
study of the transverse mass spectra, the average transverse mass (mx ) — rriQ was therefore 
calculated. To account for the small unmeasured high mx part of the kaon spectra, which 
is on a percent level only, the spectra were extrapolated by the exponential functions shown 
in Figs. [1] and |3l For the pion spectra in Pb+Pb collisions an exponential function with two 
slope parameters was used (not shown), to account for the concave shape of the distributions 
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FIG. 4. (Color online) Mean transverse mass as extracted from the transverse mass spectra at mid- 
rapidity for 7r~, K"*" and K~ versus number of participating nucleons. Open symbols show results 
from central p+p, C+C, Si+Si and Pb+Pb collisions. Results from C+C and Si+Si collisions at 
40. GeV .«.n e„.^ a. new. T.e ...... »e e....e. .o. BQ. AU o... ... poi„. 

marked by open symbols are from 3, ^, 8|. The thin vertical bars indicate the statistical errors 
which are mostly smaller than the symbols. The thick shaded bars indicate the systematic errors. 



(in logarithmic representation). An estimate of systematic errors was derived from using 
different fit ranges for the extrapolation or a single exponential for pions. The resulting 
average transverse masses and inverse slope parameters are summarized in Tables [Tll IHH 
and llVi Figure H] shows the dependence of (mx) — mo on the mean number of wounded 
nucleons (A^w) for the investigated collision systems at 40^4 GeV and 158^4 GeV beam energy. 
Results from central Pb+Pb j^, 3], p+p at 158 GeV/c 1^, C+C and Si+Si collisions js] at 
158^4 GeV are also shown. 

Owing to the mass difference and well known radial flow effect, (mx) — mo values are 
larger for kaons than for pions. For central collisions I2, 3| mean transverse masses do not 
change significantly from 40A GeV to 158A GeV. For pions, if any, only a weak centrality 
dependence is observed at 40^4 GeV and none at 158A GeV beam energy, while kaons show 
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FIG. 5. (Color online) Rapidity spectra for vr" in Pb+Pb collisions at AQA GeV and 
158^4 GeV beam energy for different centrality bins. Open symbols are reflected at mid-rapidity. 
Solid lines represent the fits of Eq. [21 The dotted and dashed lines represent HSD [36] and 



UrQMD2.3 



37l | simulations, respectively. The error bars (mostly not visible) indicate the sta- 



tistical uncertainty. 



an increase of mean transverse masses towards central Pb+Pb collisions. This increase 
is particularly pronounced when comparing to the much smaller C+C and Si+Si collision 
systems with (A^w) < 60. 



B. Rapidity spectra 

Transverse momentum spectra were measured in bins of rapidity in order to extract the 
rapidity dependence of -integrated yields. The measured px -spectra were extrapolated 
into unmeasured regions by a single exponential for both kaons and pions in C+C and Si+Si 
interactions. A double exponential had to be used for pions in Pb+Pb collisions because of 
the concave shape of their px -spectra. Rapidity distributions of px -integrated yields for the 
particles, collision systems and energies under study are shown in Figs. [5l [6l [TJ and [3 They 
are well described by two Gaussians of equal widths a which are displaced symmetrically 



15 



T3 



30 



20 



10 



20 
15 
10 
5 



■ 1 ■ ■ ■ 1 ■ 1 1 1 1 

L C4 . 

: 158AGeV 

■ HSD " 
- UrQMD 


C3. 


C2. 

9 \ 

/ \ 


. , . . . 1 . . 
/J 


C1. 

V 


CO: 

?■•/ V* 
6i % 
1 \ 

< > 


■ 40A GeV 


1 1 1 1 1 1 1 1 1 r" 


1 1 1 1 1 1 1 1 1 — |— 


— 1 — 1 1 1 1 1 1 1 

■ n \ 
1 







-2 2 -2 2 -2 2 -2 2 -2 2 

y 

FIG. 6. (Color online) Rapidity spectra of K"*" in Pb+Pb collisions at 40^ GeV and 
158^4 GeV beam energy for different centrality bins. Circles are values from the dE/dx analy- 
sis, the stars give the result from the combined dE/dx and TOF analyses. Open symbols are 
reflected at mid-rapidity. Solid lines represent the fits of Eq. [2j The dotted and dashed lines rep- 
resent HSD 36|] and UrQMD2.3 33] simulations, respectively. The error bars (mostly not visible) 



indicate the statistical uncertainty. 



around mid-rapidity by a constant Uq: 



dN 
dy 



C 



exp 



{y - yof 



+ exp 



jy + yof 

2^2 



(2) 



Using this functional form the measured rapidity spectra can be extrapolated into the un- 
measured regions, and particle yields in full phase space can be extracted. For C+C and 
Si+Si collisions at 40A GeV beam energy no combined dE/dx and TOF measurement is 
available at mid-rapidity. The resulting uncertainties in the shapes and integrals of the ra- 
pidity distributions of kaons were evaluated by applying three different fit functions (Eq. |2]): 
with two free parameters {ay and yo), with these parameters fixed to values obtained from 
fits to the corresponding distributions from central Pb+Pb coUsions at AOA GeV, and with 
a single Gaussian with its root mean square (RMS) again taken from the 40^4 GeV Pb+Pb 
data. These fits are indicated in Fig. [8] by the solid, short-dashed and dashed lines and were 
used to extrapolate into the unmeasured rapidity regions. Differences in the total yields 
between extremes are not exceeding 12% and are included in the systematic errors. The 
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quoted total and mid-rapidity yields are the means of the values obtained from the three fits. 



The widths of the rapidity distributions are quantified by RMSy = y + i/q. These values 
are tabulated together with mid-rapidity and full phase space yields in Tables [TTl IIll and 
llVi Figure [9] presents the system-size dependence of the widths of the rapidity distributions. 
While kaons show no significant change of the width as a function of the number of wounded 
nucleons, rapidity distributions of pions are significantly wider in peripheral than in central 
Pb+Pb collisions. In fact, it turns out that the width of the pion rapidity distribution 
in central Pb+Pb collisions is the same as in C+C and Si+Si reactions as well as in p+p 
interactions at 158 GeV/c. Only non-central Pb+Pb collisions deviate from this common 
behavior exhibiting a wider distribution. The widening of the pion rapidity distribution for 
peripheral collisions is probably due to pion production in interactions of participants with 
spectator matter which had been already invoked to explain the proton rapidity spectra in 
the same data sets jlf. 



>,20 
T3 



T3 15 
10 
5 



. 1 . . . 1 . . . , . 

C4 

■ 158AGeV " 

: ■ HSD 
- UrQMD 


C3 

< 1 1 1 1 1 1 1 1 1 > 


C2 ; 

i \ ' 

-1 1 1 1 1 1 1 1 1 1 > 


1 1 1 1 ■ 1 ■ ■ ■ 1 ■ 

C1 

Si \ 
4' ^ 

i \ 

/ \ 
^1 1 1 1 1 1 1 1 1 1^ 


r \ 

-7 \ - 
I' \ ■ 

1 1 1 1 1 1 1 1 1 


40A GeV 











10 



FIG. 7. (Color online) Rapidity spectra of K~ in Pb+Pb collisions at 40A GeV and 158A GeV beam 
energy for different centrality selections. Circles are values from the dE/dx analysis, the star gives 
the result from the combined dE/dx and TOF analysis. Open symbols are reflected at mid' 



36| 



rapidity. Solid lines represent the fits of Eq. [2j The dotted and dashed lines represent HSD 
and UrQMD2.3 {s?] simulations, respectively. The error bars (mostly not visible) indicate the 
statistical uncertainty.. 
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FIG. 8. (Color online) Rapidity spectra for vr^ and for semi-central C+C and Si+Si collisions at 
40^4 GeV beam energy. Solid lines show free fits of double Gaussians. Open symbols are reflected 
at mid-rapidity. For also double Gaussians with width a and shift yo as for central Pb+Pb 
collisions at 40A GeV beam energy were fitted (short dashed line) as well as a single Gaussian with 
the same RMSy (dashed line). The error bars indicate the statistical uncertainty. 

C. Particle yields 



Fig. [To] shows the system-size dependence of particle yields normalized to the mean num- 
ber of wounded nucleons. Earlier data on central C+C, Si+Si and Pb+Pb collisions are 
2, 



plotted as well 



, |3|, |8[. Data for p+p collisions are taken from NA49 where available jo], [l^ 
and from parametrizations by Rossi et al. [38:|. At 40 GeV/c we include results from a new 
study of the energy dependence of kaon production reported in [lol]. Also shown are calcu- 
lations by the HSD [36\ and UrQMD2.3 37| transport models, as well as by the core-corona 
model (the latter for the Pb+Pb collision system only, see also Fig. [T2|) 39| . Normalized 
pion yields are rather independent of centrality or system size while the normalized yields of 
kaons show a steep increase for small system sizes followed by a slow rise or even saturation 
for higher centralities. The mean number of wounded nucleons {{N^)) introduces a model 
dependence in these ratios, in particular for small systems, as for the semi-central C+C and 
Si+Si collisions the calculations have large uncertainties. Yields are therefore alternatively 
normalized to the number of pions. In Fig. [11] these ratios are presented. Normalized kaon 
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FIG. 9. (Color online) Width of the rapidity distribution described by the RMS values for vr^, 
and K~ versus mean number of participating nucleons. Open symbols show results from 
p+p, C+C, Si+Si and central Pb+Pb collisions. Except for the^new results for C+C and Si+Si 
collisions at AO A GeV beam energy the data are taken from 2, 



a, Isl . Filled symbols represent the 



37t | simulations 



centrality selected Pb+Pb data. Lines show results from HSD [36] and UrQMD2.3 
as indicated in the figure. The thin vertical bars indicate the statistical errors, which are mostly 
not visible. The thick shaded bars indicate the systematic errors. 



yields at both energies exhibit the same dependence as seen in Fig. [TDl A closer look reveals 
that at 158^4 GeV the K/vr ratios in the small collision systems C+C and Si+Si seem to fol- 
low the trend set by the centrality dependence of the Pb+Pb data, whereas at AOA GeV the 
normalized kaon abundance in peripheral Pb+Pb collisions is lower than expected from the 
extrapolation of the K/tt ratio in C+C and Si+Si collisions. 
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FIG. 10. (Color online) Particle yields divided by the mean number of wounded nucleons as a 



function of centrality. Closed symjo^ 



s show results from centrality selected Pb+Pb collisions, 
open circles p+p data from NA49 [9|, [lOl] as well as C+C, Si+Si (at 158^ GeV [8]) and central 
Pb+Pb collisions from 3) ) • The stars denote the values obtained for p+p interactions from the 
parametrization "B" in reference [ss]- At 40 GeV/c they are compared to values extracted from 



371 ] calculations as 



^). The thin 



Fig. 132 reported in [10]. Lines show results from HSD [36] and UrQMD2.3 
well as from the core-corona model (the latter for the Pb+Pb collision system only, 
vertical bars indicate the statistical errors, which are mostly not visible. The thick shaded bars 
indicate the systematic errors. 
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FIG. 11. (Color online) Kaon yields divided by mean number of pions (vr^) = 0.5 • ((tt ) + (tt^)) 
as a function of centrality. Closed symbols show results from centrality selected Pb+Pb collisions, 



open symbols p+p data from NA49 



Pb+Pb collisions from 
in reference 



as well as C+C, Si+Si (at 158^ GeV fSj]) and central 



). The stars denote the values obtained from the parametrization "B" 



and UrQMD2.3 

n. 



simulations and the core- 



381 . Lines show results from HSD 
corona model (the latter for the Pb+Pb collision system only, [39]). The thin vertical bars indicate 
the statistical errors, which are mostly not visible. The thick shaded bars indicate the systematic 
errors. 



V. DISCUSSION 



A. Comparison to transport models (HSD and UrQMD) 

The centrality of the events generated by the models was derived from the corresponding 
fraction of the total (geometrical) cross section with the impact parameter as order parame- 
ter. The corresponding (A^„) values were those derived from the experimental data for HSD. 
For the UrQMD events the number of nucleons in the nuclear overlap volume (A^„) was cal- 
culated event-by-event based on the chosen impact parameter. Pion rapidity distributions 
(Fig. E]) are well reproduced by both models at iOA GeV, whereas their predictions give 
slightly higher yields than the data at 158A GeV. For K"*" mesons (Fig. |6]) the distributions 
from both models are 20% low at 40^4 GeV. For kaons at 158yl GeV HSD calculations agree 
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well with experiment, whereas UrQMD2.3 underestimates our data. The overall agreement 
betweeen model calculations and data for the K~ mesons (Fig. [7]) is generally good for 
both models, even though there is a small underestimation of the data by UrQMD2.3 at 
158^4 GeV and by HSD at 40^4 GeV. No model data were provided for comparison to our 
C+C and Si+Si data. As to the shapes of the rapidity distributions both models describe 
well the increase of RMSy with increasing impact parameter at both energies (Fig. [9]). This 
increase is most (least) pronounced for the pions (K~). We interpret this rise as due to 
interactions of fireball particles with spectator remnants. 

Normalized multiplicities are studied in Figs.[TO]and[TTl (A^) / {N^) for pions is reproduced 
by UrQMD2.3 at both energies. HSD is 10% high, which is due to larger widths rather than 
higher amplitudes of the corresponding dN/dy distributions (cf. Fig. E]). We now turn to 
the normalized yields of the K"*" meson. At 40^4 GeV the results from both models are 
low by 20%. At 158A GeV HSD results agree with our data, wheras UrQMD2.3 results 
underpredict them. The K~ yields are well desribed by both models at 40^4 GeV. At 
158A GeV K~ mesons yields are reproduced by HSD, whereas they are underpredicted by 
UrQMD2.3 as was the case for K"*" yields. The calculations of both models for the kaon to 
pion ratios yield values which are lower than the experimental data at both energies (see 
Fig. [TT]) . although for different reasons. While UrQMD2.3 mainly underpredicts the kaon 
yields, HSD overpredicts the pions yields (see Figs. [TOj) . 



B. Transverse mass distributions and thermal freeze-out 

Both the thermal freeze-out temperature and radial flow affect the observed mx distri- 
butions of produced particles. The centrality dependence of (mx) —rriQ (near mid-rapidity), 
which is shown in Fig. HI exhibits a clear trend of an increase of (mx) — mo with centrality 
(for Pb+Pb) and system size (p+p and light versus heavy systems), which is most pro- 
nounced for kaons of both charges at 158A GeV. This finding is consistent with the picture 
that nuclear stopping and thus energy deposition increases with increasing centrality and 
from small to large collision systems, giving rise to higher temperatures or collective radial 
flow or both. 
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C. System size dependence of the K/vr ratio 
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FIG. 12. (Color online) Core fraction obtained in the core-corona model for different collision 
systems as function of the number of participating nucleons [39]. The dashed line connects the 
endpoints of the curves for the different collision systems. /((A^w)) denotes the fraction of those 
nucleons, which undergo more than one collision. 

The system size dependence of the K/vr ratios (see Fig. [TT]) hints at a behavior expected 



for the core-corona ansatz 39|. However at 40^4 GeV there are significant deviations of the 
core-corona calculations from the experimental data. In peripheral Pb+Pb collisions the 
data tend to be below the model values, most pronounced for K~, and the ratios measured 
for the small collision system (C+C and Si+Si) tend to be above the values calculated for 
the Pb+Pb system. The latter finding is in line with the model prediction for different 
system sizes as shown in Fig. [121 H shows the core fraction obtained in the core-corona 
model for different collisions systems as function of the number of participating nucleons 



39| . Central collisions of small projectiles reach their maximum already for smaller volumes 
(number of wounded nucleons) than heavy systems like Pb+Pb. At 158^4 GeV (see right 
column of Fig. [TTj) the data points from peripheral Pb+Pb collisions follow nicely the core- 
corona prediction for the Pb+Pb system. Also the C+C and Si+Si points are consistent 
with this prediction (for Pb+Pb), although according to Fig. [12] they should be higher even 
for non-central collisions of the small C+C and Si+Si systems. However, here the core- 
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corona approach is applied to small system sizes for which the simplification of separating 
the collision zone into one central cluster and p+p collisions in the periphery might not 
be valid anymore. Instead, percolation calculations suggest that for such collision systems 



several clusters of smaller size will be formed 
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FIG. 13. (Color online) Centrality dependence of the double ratio of the K"^ to vr^ multiplicities 
scaled to the ratios from the corresponding most central measurement minimum bias Pb+Pb 
lisions at SIS 11] and AGS [3) [3] energies (a), at SPS energies (b) and at RHIC 



(Au+Au) co: 

n 

energies |14l . 



isl (c). The dashed lines in panels (a) and (c) are meant to guide the eye and are 
also drawn as a reference in panel (b). Only statistical errors are shown. 



In Fig. US] we summarize the K~^/tt^ ratios for centrality selected Au+Au (Pb+Pb) col- 
lisions at different beam energies. The ratios at the different centralities were scaled to the 
ratios measured in the most central collisions. The observed saturation of strangeness pro- 
duction at SPS energy and above can be understood in the context of statistical models by 
the approach to a grand-canonical description. Increasing the energy available for particle 
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production does no longer change the multiphcities of strangeness carrying particles relative 
to pions 3|, which means that strangeness becomes fully saturated. Saturation already 
sets in for rather peripheral collisions. The reduction of relative strangeness production 
towards even smaller system sizes has been discussed in the framework of the core-corona 
model in the previous paragraph. An implementation of this ansatz into the statistical 
model and a detailed comparison of model calculations to RHIC data can be found in 25 |. 
In order to describe the data satisfactorily it is essential to use a volume dependence which 
is not just proportional to (A^„). The situation is different at low (SIS) energies. Here the 
K^/tt^ ratio does not saturate and its dependence on the number of wounded nucleons is 
rather linear. In contrast to the situation at high energies, statistical models can describe 
the observed centrality dependence at low energy within a canonical ensemble and with a 
volume proportional to the number of wounded nucleons V oc (A^w) |20|. Such a behavior 
would also be expected in a hadronic rescattering scenario as realized in microscopic models 
13| . The results obtained at 40 A GeV are intermediate between those at top SPS and AGS 
energies. We thus observe a smooth transition between both scenarios which occurs in the 
SPS energy region. 



D. Es as function of system size 

The strangeness to entropy ratio or total relative strangeness production described by the 
Wroblewski factor = 2{ss) / {uu + dd) can be approximated by the measurable quantity 
Es= {Na + 2[Nk+ + A^i^-])/(1.5[A'^+ + A^^-]) |2|]. The data presented in this paper com- 
bined with earlier ones published in jol, IsJ allow to study the dependence of Es on system 
size for 40A GeV and 158A GeV beam energy with the 158A GeV data covering the whole 
range of possible sizes. This dependence is shown in Fig. [H] and can now be compared 
directly to the dependence obtained from statistical model calculations assuming for sim- 
plicity that Es (X v, with rj being the canonical strangeness suppression factor, see for 



example Refs. [l8|, ll9|. In statistical models the relative strangeness production depends 



on the volume of the ensemble. Assuming naively that the volume is proportional to the 



number of wounded nucleons ( 18|, ll9|) yields an increase of relative strangeness production 
(dashed line in Fig. [T^ which is significantly steeper than the increase observed in the data. 
On the other hand, the smooth rise of Es with system size and the saturation around 60-100 
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FIG. 14. (Color online) Total relative strangeness approximated by Eg for centrality selected 
minimum bias Pb+Pb collisions at 40A GeV (squares) and 158A GeV (circles) as function of 



40|] and central Pb+Pb 



(Nw)- The open symbols are from p+p [9|, [lOj, C+C, Si+Si ^l, S+S 
4] at 158j4 GeV. The filled symbols represent the data of centrality selected Pb+Pb collisions at 
AO A GeV and 158^ GeV (this analysis and {3]). The lines indicate calculations within the statistical 
model assuming proportionality between ensemble volume and (A^w) (dashed), with volumes as 
derived from percolation calculations for small systems (dotted line) and Pb+Pb collisions (solid 
line) [21], as well as calculations based on the core-corona ansatz at 158^4 GeV (dash-dotted line). 
Only statistical errors are shown. 



wounded nucleons is reproduced within experimental errors by a percolation model [2l|]. In 
this model the volumes relevant for strangeness production by statistical hadronization are 
calculated within a percolation ansatz. This yields several smaller clusters in collisions of 
light nuclei or in peripheral Pb+Pb collisions, and a large cluster accompanied by small 
clusters (p+p interaction like) in the periphery for central Pb+Pb collisions. The latter 
finding justifies the simplified ansatz used in the core-corona approach, which is shown by 
the dash-dotted line. The functional dependences change slightly between light and heavy 
nuclei (dotted and solid lines in Fig. [T^ . both describing the data at 158 A GeV beam energy 
well. The Pb+Pb data show a weaker increase at 40^4 GeV than at 158^4 GeV. 
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VI. SUMMARY 



The NA49 collaboration measured the rapidity and transverse mass distributions of kaons 
and pions in Pb+Pb coUisions at different centralities. These data are complemented by 
results from the small systems C+C and Si+Si. We find that the centrality dependence of 
kaon and pion production is reproduced by microscopic transport model calculations (HSD 
and UrQMD2.3) within 20% or better. The mean transverse mass evolves with centrality 
as expected for increased stopping and energy deposition. The centrality dependence of the 
widths of the rapidity distributions in Pb+Pb collisions does not smoothly connect to results 
from the small systems p+p, C+C, and Si+Si. We attribute this behavior to reinteractions 
of produced particles with spectator remnants, which is possible only in the heavy collision 
system. In Pb+Pb coUisions the K/tt ratios show a smooth increase with centrality at both 
beam energies with saturation setting in around 100-200 (60-100) wounded nucleons for 
40 A GeV (158A GeV). In particular at 40 A GeV beam energy the ratios measured in C+C 
and Si+Si collisions tend to be higher than the ratios in peripheral Pb+Pb interactions. The 
strong rise for small systems can be understood in the framework of statistical models, if the 
relevant volume is not proportional to the number of participants but is a superposition of 
relatively small subvolumes for small system sizes and one large volume in central collisions 
accompanied by p+p like interactions in the corona. 
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1 


22.72 ±0.31 ±2.3 


78.5 ±2.8 ±7.8 


1.17 ±0.02 ±0.06 


302 ±5 ±15 


229 ±3 ±12 




2 


15.86 ±0.23 ±1.6 


57.1 ±2.9 ±5.7 


1.18 ±0.02 ±0.06 


300 ±5 ±15 


227 ±2 ±12 




3 


10.01 ±0.17 ±1.0 


34.9 ±2 ±3.5 


1.21 ±0.02 ±0.06 


283 ±6 ±14 


215 ±3 ±11 




4 


6.39 ±0.13 ±0.6 


23.2 ±1.9 ±2.3 


1.22 ±0.02 ±0.06 


277 ±7 ±14 


209 ±3 ±11 


K- 





16.8 ±0.43 ±1.7 


53.95 ±2.15 ±4.3 


1.12 ±0.03 ±0.06 


298 ±9 ±15 


225 ±5 ±11 




1 


13.38 ±0.24 ±1.3 


43.18 ±2.03 ±3.5 


1.15 ±0.03 ±0.06 


294 ±7 ±15 


221 ±4 ±11 




2 


9.54 ±0.17 ±0.9 


31.32 ±2.03 ±2.5 


1.25 ±0.1 ±0.06 


285 ±8 ±14 


217 ±3 ±11 




3 


6.26 ±0.12 ±0.6 


20.45 ±1.69 ±1.6 


1.29 ±0.03 ±0.07 


273 ±8 ±14 


210 ±4 ±11 




4 


3.92 ±0.11 ±0.4 


12.24 ±1.18 ±1.2 


1.21 ±0.07 ±0.06 


264 ±10 ±1.1 


202 ±5 ±10 



TABLE II. Summary of data for pion and kaon production in centrality selected Pb+Pb collisions at 
158A GeV beam energy: rapidity density dN/dy at mid-rapidity (tt": < y < 0.4, K='=: average of 
result from TOF analysis (—0.2 < y < 0.2) and dE/dx analysis (0.1 < y < 0.3)), total multiplicity 
(A^), RMS width of the rapidity distribution RMSy, average transverse mass at mid-rapidity 
(mx) — mo, inverse slope parameter of mx spectra at mid-rapidity (vr": < y < 0.2, fitrange 0.2 
< mx - mo < 0.7 GeV/c^; K=^: -0.2 < y < 0.2 (TOF analysis), fit range < mx - mo < 1.0 
GeV/c^). The first error is statistical, the second systematic. 
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Centr. 


diV/dy 


(iV) 


RMSy 


(mT) — rriQ 


T 




class 








MeV/c2 


MeV 


TT 





118 ±0.6 ±12 


348 ±2.4 ±17.4 


1.11 ±0.01 ±0.06 


270 ±3 ±27 


172 ±3 ±9 




1 


93.8 ±0.5 ±9.4 


286 ±2.4 ±14.3 


1.15 ±0.01 ± 0.06 


272 ±3 ±27 


172 ±3 ±9 




2 


67.8 ±0.4 ±6.8 


214 ±2.4 ±10.7 


1.19 ±0.01 ± 0.06 


273 ±2 ±27 


178 ±2 ±9 




3 


43.7 ±0.3 ±4.4 


146.3 ±2.2 ±7.3 


1.25 ±0.01 ±0.06 


269 ±3 ±27 


172 ±3 ±9 




4 


29.1 ±0.2 ±2.9 


101.1 ±1.9 ±10 


1.29 ±0.01 ± 0.06 


263±3 ±26 


168 ±3 ±9 


K+ 





20.1 ±0.6 ±1.9 


60.0 ±4.7 ±4.8 


1.00 ±0.02 ±0.05 


312 ±13 ±15 


238 ±7 ±12 




1 


16.2 ±0.4 ±1.5 


49.1 ±4.3 ±3.9 


1.05 ±0.02 ±0.05 


307 ±11 ±15 


238 ±6 ±12 




2 


10.9 ±0.3 ±1.0 


34.2 ±3.5 ±2.7 


1.03 ±0.01 ±0.05 


309 ±12 ±15 


234 ±6 ±12 




3 


6.9 ±0.2 ±0.7 


22.0 ±2.8 ±1.7 


1.06 ±0.02 ±0.05 


291 ±14 ±15 


221 ±7 ±11 




4 


4.2 ±0.2 ±0.4 


13.6 ±2.3 ±1.4 


1.08 ±0.03 ±0.05 


265 ±16 ±13 


208 ±9 ±11 


K- 





8.5 ±0.3±0.7 


21.0 ±3.1 ±2 


0.89 ±0.03 ±0.05 


323 ±19 ±15 


239 ±10 ±12 




1 


6.2 ± 0.2±0.5 


15.9 ±2.7 ±1.3 


0.91 ±0.03 ±0.05 


294 ±16 ± 15 


231 ±9 ±12 




2 


4.7 ±0.2±0.4 


11.4 ±2.3 ±0.9 


0.89 ±0.04 ±0.05 


304 ±16 ±15 


227 ±8 ±12 




3 


2.7 ±0.1 ±0.2 


7.0 ±1.8 ±0.6 


0.88 ±0.03 ±0.05 


275 ±20 ±14 


207 ±10 ±11 




4 


2.0 ±0.1±0.2 


4.6 ±1.5 ±0.5 


0.85 ±0.04 ±0.05 


291 ±24 ±15 


219 ±14 ±11 



TABLE III. Summary of data for pion and kaon production in centrality selected Pb+Pb colli- 
sions at 4QA GeV beam energy: rapidity density dN/dy at mid-rapidity {tt": < y < 0.4, K='=: 
result from TOP analysis (-0.1 < y < 0.1)), total multiplicity (iV), RMS width of the rapidity 
distribution RMSy, average transverse mass at mid-rapidity (mx) — mo, inverse slope parameter 
of niT spectra at mid-rapidity (tt": < y < 0.2, fitrange 0.2 GeV/c^ < rriT — mo < 0.7 GeV/c^; 
K=^: -0.1 < y < 0.1 (TOP analysis), fit range < mx - mo < 0.8 GeV/c^). The first error is 
statistical, the second systematic. 
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Collision 


dN/dy 


(N) 


RMSy 
y 


(mT) — mn 
MeV/c^ 


T 
MeV 


TT 


C+C 


3.25 ±0.06 


10.0 ±1.4 


1.18 ±0.04 ±0.1 


247 ±11 ±25 


169 ±8 ±9 




Si+Si 


11.40 ±0.19 


34.6 ±2.7 


1.16 ±0.03 ±0.1 


243 ±10 ±24 


162 ±6 ±9 


K+ 


C+C 


0.46 ±0.06 


1.26 ±0.78 


1.10 ±0.11 ±0.15 








Si+Si 


1.82 ±0.30 


5.1 ±1.2 


1.01 ±0.15 ±0.15 






K- 


C+C 


0.19 ±0.02 


0.44 ±0.4 


0.80 ±0.08±0.12 








Si+Si 


0.87 ±0.07 


2.0 ±1.0 


0.89 ±0.09±0.12 







TABLE IV. Summary of data for pion and kaon production in C+C and Si+Si collisions at 
40^ GeVbeam energy: rapidity density dN/dy at mid-rapidity (vr": < y < 0.4, K^: y = 
average of the three fit functions), their multiplicities {N), and the RMS widths RMSy of their 
rapidity distributions. For pions also average transverse masses (mT) — rriQ and the inverse slope 
parameters T of the mx spectra are given (0 < y < 0.2, fitrange 0.2 < mx — mo < 0.7 GeV/c^). 
The first error is the statistical, the second systematic. The systematic errors on the yields are 
15%. 
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